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ABSTRACT 

We discuss two mistreatments of damped Lya (DLA) kinematic analysis that were first performed 
by Haehnelt, Steinmetz, & Rauch (1998; hereafter HSR98) and have recently been repeated by Hong et 
al. (2010; astro-ph/1008.4242_vl,v2; hereafter H10). Each mistreatment led to the improper excising 
of simulated absorption profiles. Specifically, their analyses are strictly biased against DLA sightlines 
that have low HI column density Nm < 10 20 ' 5 cm -2 , very high Nm values, and (for all iVjji) sightlines 
with low Av (< SOkms" 1 for HSR98; < 20 — 30kms _1 for H10). None of these biases exist in 
the observational analysis. We suspect these mistreatments compromise the results that followed. 
Hopefully this posting will prevent their repetition in the future. 

Subject headings: galaxies: evolution — intergalactic medium — quasars: absorption lines 



1. INTRODUCTION 

In 1997, we first measured the kinematic characteris- 
tics of a small sample (N — 17) of damped Lya sys- 
tems (DLAs) as a mea ns to assess the dynamica l mo- 
tions of high z galaxies (|Proc haska & Wolfe 1997b, here- 
after PW97). We compared these data against simplistic 
models of protogalaxies (e.g. rotating disks, infalling gas 
clumps) and favored a thick (h > 1 kpc) disk scenario 
with rotation speeds v c ~ 200 km s _1 . Furthermore, we 
ruled out standard CD M scenarios that c onsidered DLAs 
as rotating disks (e.g. IKauffmannlll996f ). These conclu- 
sions hold today. 

The following year, IHaehnelt et al.l ([1998., hereafter 
HSR98) published a competing model within the CDM 
framework that described DLAs a merging protogalactic 
clumps. Their hybrid analysis combined results from a 
small set of 'zoom-in' numerical simulations with a Press- 
Schecter treatment of the dark matter (DM) halo mass 
function, and a heuristic assumption for the DM halo 
dependence of DLA gas cross-section. HSR98 concluded 
that this model could reproduce the DLA kinematics pro- 
vided that low mass DM halos (v c K SOkms" 1 ) [Check; 
could be 50] did not host DLAs (sec Barnes fc Haehnelt] 
2009, for a modern update of this model). 

Many other models have since been proposed with a 
range of success ( e.g. iNulsen et all 119981 : iMaller et al 



2001b lSchavell200l1: iRazoumov et al.ll2 008: Po ntzen et ah 
20081) . Most recently, |Hmig_et_aL| ^QM hereafter H10) 



proposed that galactic-scale winds are an essential ingre- 
dient to explaining DLA kinematics. In reviewing this 
paper, we noted two treatments of the simulated spectra 
that depart from the prescription applied to the DLA ob- 
servations (as detailed in PW97). To our surprise, we dis- 
covered that these mistreatments were first introduced by 
HSR98. We expect that H10 simply repeated the HSR98 
analysis without realizing the errors. As H10 themselves 
demonstrate, these have important consequences for the 
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Figure 1. Low-ion transitions for th e DLA at z = 1.920 toward 
Q2206-19 UProchaska fc Wolfe! ri997aj ). Note how closely the var- 
ious transitions trace one another, whether or not they arise from 
the same ion. Analysis of any of the lines that satisfy the criterion 
in Equation Q] (only Si II 1808, Ni II 1741 do) will yield nearly iden- 
tical kinematic characteristics. At this S/N, one could even include 
weaker transitions (e.g. Ni II 1751, Cr II 2056). We have never pre- 
cluded kinematic analysis of a DLA based on this criterion; it only 
serves to determine which transition is used. 
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results derived from the simulations. With this posting, 
we wish to highlight these issues in the hopes of prevent- 
ing their further propagation. 

2. THE DEPTH CRITERION 
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Because the H I Lyman series of DLAs is highly satu- 
rated, we cannot assess the gas kinematics with neutral 
hydrogen. Instead, DLA kinematics have been assessed 
using low-ion, metal-line transitions (e.g. Si II 1808). 
To minimize the effects of line-saturation and statisti- 
cal noise in the spectra, we imposed a criterion on the 
depth of the transition that is analyzed: 

0.1 < min (j-^) - °' 6 W 

where I is the observed intensity and I c is the continuum 
of the quasar local to the transition. This criterion pre- 
cludes lines that are heavily saturated (i.e. I/I c < 0.1) 
also weak profiles to focus the analysis on the bulk of 
the gas while maintaing a good S/N ratio (see PW97 for 
more details) . Conveniently, nature provides enough UV 
transitions at A < 2000A such that we have always iden- 
tified at least one low-ion transition for analysis in every 
DLA observed, no matter its H I column density and 
metallicity. Furthermore, the low-ion transitions very 
precisely trace one another (Figure □ Procha]S3[200i 
see also) such that all lines satisfying this criterion have 
nearly identical kinematic characteristics. To emphasize, 
no DLA has been precluded from kinematic analysis be- 
cause of this criterion. Indeed, the original PW97 anal- 
ysis relied on transitions with a wide range of Xf values 
(e.g. Fe II 1608, Sill 1808). 

HSR98 generated metal-line absorption profiles from 
their numerical simulations (i) by assuming the H I 
gas had a metallicity [Si/H] =—1 everywhere and (ii) 
by adopting the atomic data for the Si II 1808 transi- 
tion. These authors then only considered profiles "sat- 
isfying the criterion in equation ([lj". In other words, 
the Si II 1808 profiles that violated this criterion were 
excised from the analysis. In fact, this likely includes 
nearly all of the profiles that they presented in their Fig- 
ures 3-6 as representative of their analysis. H10 followed 
the same steps as HSR98 and thereby excised 30 — 60% 
of the profiles from analysis. 

It is impossible for us to precisely predict the implica- 
tions of this mistreatment, but there are some obvious 
biases. Figure [2] shows the Si II 1808 profile for a DLA 
at the H I column density threshold TVhi = 10 2a3 cm~ 2 
with metallicity [Si/H] = —1, and a single component 
with Doppler parameter b = 8kms _1 . We also indi- 
cate the depth criterion expressed in Equation [TJ In a 
mistreatment of this criterion, the sightline would have 
been excised from kinematic analysis. In a proper treat- 
ment, one would identify another, strong transition (e.g. 
Si II 1526, Fe II 1608) to perform the analysis. In short, 
the mistreatment biases analysis against sightlines with 
low Nhi value. Figure [2] also shows that sightlines with 
very large -/Vhi value may produce saturated profiles that 
also violate the depth criterion. Therefore, this mistreat- 
ment also biases the results against sightlines with very 
high Nja. Of course, nature provides many more low ./Vhi 
sightlines than high iVm~ sightlines such that the overall 
bias is against the former. The effects kinematically are 
more difficult to predict. 

3. A Av THRESHOLD 

The other mistreatment of the simulated profiles by 
HSR98 was: "For the discussion of velocity widths a 
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Figure 2. Simulated Si II 1808 profiles for (blue) a DLA with 
./Vhi = 10 20 ' 3 cm -2 , a metallicity [Si/H] = —1, and a single com- 
ponent with Doppler parameter b = 8 km S ; and (red) LA with 
total JVhi = 10 21 ' 3 cm -2 , a metallicity [Si/H] = —1, and two com- 
ponents each with b = 8kms _1 . The shaded region shows the 
depth criterion described by Equation [T] Neither of these profiles 
satisfy the criterion: the low Nm system is too weak and the high 
Mm has too saturated of a Si II 1808 profile. These profiles would 
have been excised in the analysis of HSR98 and H10. In the obser- 
vational analysis, however, we would have identified other low-ion 
transitions that do satisfy the criterion and would not have excised 
the sightline from kinematic analysis. 

minimum threshold of Av > 30kms _1 was imposed 
on both the observed and simulated velocity widths to 
avoid incompleteness effects." . Observationally, one has 
never imposed such a cut. There is no equivalent to the 
magnitude limit in absorption line studies, especially for 
kinematic analysifl Every DLA known that has been 
observed at even moderate S/N and spectral resolution 
shows several low-ion, metal transitions. Indeed, this has 
established the so -called 'floor' in metalli c ity at ~ 1/1000 
solar for DLAs (jProchaska et al.l 120031 : IPcnpras e et al.l 

HoTJ). 

Imposing a Av < 30 km s _1 threshold actually has a 
minor effect on the DLA observations because so few sys- 
tems exhibit such low Av values (see iWolfe et al.l I2005L 
for the latest compilation). In contrast, models may pro- 
duce a preponderance of systems with small Av values. 
Adopting a Aw threshold, therefore, can greatly affect the 
results when one compares the model predictions against 
the observations. Indeed, the KS probability of the fa- 
vored wind model in H10 decreases from 39% (fine agree- 
ment) to 2.8% (nearly ruled out) as they reduced the 
threshold from 30 km s _1 to 20 km s . It is quite pos- 
sible that their model would be ruled out at > 99% c.l. 
if no threshold were imposed (i.e. the proper treatment). 
We can only speculate on the implications for HSR98. It 
seems likely that the primary result from their analysis 
of the numerical simulations, that Av ~ 0.6w c , would be 
greatly diminished (e.g. examine their Figure 10a). Pre- 
sumably, this could be tested in the 'no wind' models of 

3 It has been oft debated that a limit exists for metallicity mea- 
surements owing to dust obscuration by the DLA of the background 
quasar. This would almost certainly bias one against sightlines 
with large velocity width, not small. 
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H10. 

4. CONCLUDING REMARKS 

Before concluding, we (sheepishly) apologize for not 
having identified these issues previously. We also fully 
acknowledge that the many successes of CDM will not 
be overturned by simple mistreatments of DLA kine- 
matic analysis. Nevertheless, the results appear to re- 
main a challenge to modern theories of galaxy formation. 
Hopefully the tension that lies between the observations 
and theoretical prediction will reveal key aspects of ISM 
physics and feedback in the early universe. 

JXP and AMW are supported by NSF grant (AST- 
0709235). JXP acknowledges the hospitality of UC San 
Diego and CASS where he is spending a sabbatical. 
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